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Abstract. We consider a class of translationally invariant magnetic fields such that 
the corresponding potential has a constant direction. Our goal is to study basic spec- 
tral properties of the Schrodinger operator H with such a potential. In particular, 
we show that the spectrum of H is absolutely continuous and we find its location. 
Then we study the long-time behaviour of solutions exp(— iHi)/ of the time depen- 
dent Schrodinger equation. It turnes out that a quantum particle remains localized in 
the plane orthogonal to the direction of the potential. Its propagation in this direction 
is determined by group velocities. It is to a some extent similar to a evolution of a 
one-dimensional free particle but "exits" to +00 and —00 might be essentially different. 



1. Introduction 

1.1. Traslationally invariant magnetic fields B{x) = b2{x), b^^x)), x = {xi,X2, X3) 

div B{x) = 0, give important examples where a non-trivial information can be obtained 
about spectral properties of the corresponding Schrodinger operators H. We suppose 
for definiteness that B{x) does not depend on the xs-variable so that H commute with 
translations along the xa-axis. There are two essentially different (and in some sense 
extreme) classes of traslationally invariant magnetic fields. 

The first class consists of fields B{x) = (0, 0, 63(xi, X2)) of constant direction. For 
such fields, the momentum p of a classical particle in the xs-direction is conserved, 
and in the Schrodinger equation the variable X3 can be separated. Thus, we arrive to 
a two-dimensional problem in the (xi, X2)-plane. Furthermore, if 63 is a function of 
r = {xl + x'^Y^'^ only, then we get a set of problems on the half-line r > labelled 
by the magnetic quantum number m. The most important example of this type is a 
constant magnetic field 63 (r) = const (see [ID])- Some class of functions 63 (r) decaying 
as r — > 00 was discussed in P2] (see also [2]) where new interesting effects were found. 
Another famous case is &3(xi,a;2) = 6{xi,X2) (5(-) is the Dirac delta-function) studied 
in [1]. Scattering by an arbitrary short-range (decaying faster than e > 0, as 

|x| — >• 00) magnetic field 63(xi,X2) turns out to be rather similar to this particular case 
(see [16j). 

The second class consists of fields B{x) = (&i(a;i, X2), &2(a^i, 2:2), 0) orthogonal to the 
X3-axis. In this case the corresponding magnetic potential A{x), defined (up to gauge 
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transformations) by the equation cmlA{x) = B{x), can be chosen as 

= (0,0,-a(xi,X2)) (1.1) 

so that it has the constant direction. In contrast to fields of the first class, now the 
variable X3 cannot be separated in the Schrodinger equation. Nevertheless due to the 
invariance with respect to translations along the Xs-axis the operator (we always suppose 
that the charge of a particle is equal to 1) 

H={iV + A{x)Y (1.2) 

can be realized, after the Fourier transform in the variable X3, in the space L^(M; L^(M^)) 
as the operator of multiplication by the operator-valued function 

H{p) = -A+{aixi,X2)+py:L\R'')^L\R''), peR. (1.3) 

Moreover, if a{xi,X2) = a{r), then the subspaces with fixed magnetic quantum num- 
ber m G Z are invariant subspaces of H{p) so that the operator H{p) reduces to the 
orthogonal sum over m G Z of the operators 

1 (i d Tfi^ 

Hrnip) = ~-^r— + — + (a(r) + p^ (1.4) 
r dr dr 

acting in the space Ti = L^(M+; rdr). In this case the field is given by the equation 

5(x) = 6(r)(-sin^,cos^,0) (1.5) 

where b{r) = a'(r) and 6 is the polar angle. Thus, vectors B{x) are tangent to circles 
centered at the origin. An important example of such type is a field created by a current 
along an infinite straight wire (coinciding with the xs-axis). In this case b{r) = bQr~^ so 
that a(r) = bQlnr. The Schrodinger operator with such magnetic potential was studied 
in [15j. 

1.2. In this article we consider magnetic fields fll.Sp with a sufficiently arbitrary func- 
tion b{r). Our goal is to study basic spectral properties of the corresponding Schrodinger 
operator H such as the absolute continuity, location and multiplicity of the spectrum, 
as well as the long-time behaviour of the unitary group exp (— zHt). We emphasize that 
for magnetic fields considered here, the problem is genuinely three-dimensional, and 
actually the motion of a particle in the Xs-direction is of a particular interest. 

Using the cylindrical invariance of field (11.51) . we can start either from translational 
or from rotational (around the Xs-axis) symmetries. The rotational invariance implies 
that the operator H is the orthogonal sum of its restrictions on the subspaces of 
functions with magnetic quantum number m G Z = {0, ±1, ±2, . . .}. It can be identified 
with the operator (we keep the same notation for this operator) 

l^m = --^r^ + ^ + {^-^ + a{r)r (1.6) 
r dr dr 0x3 

acting in the space = L^(M+ x M; rdrdxs). In view of the translation invariance, every 
operator can be realized (again after the Fourier transform in the variable X3) in 
the space L^(]R; L^(M+; rrfr)) as the operator of multiplication by the operator- valued 
function Hm{p) defined by (11.41) . 
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Suppose now that 6(r) does not tend to zero too fast so that a(r) ^ oo as r — > oo. 
Then the spectrum of each operator Hmip) is discrete. Let A„^m(p), ?7,gN = {1,2, ...}, 
be the increasing sequence of its eigenvalues (they are simple and positive), and let 
4'n,m{r,p) be the corresponding sequence of its eigenfunctions. The functions Xn,m{p) 
are known as dispersion curves of the problem. They determine the spectral properties 
of the operator Hm- 

Note that if a(r) is replaced by —a{r), then \n,m{p) is replaced by —\n,m{p), so that 
the case a(r) — oo as r ^ oo is automatically included in our considerations. Recall 
that, for a magnetic field B{x), the magnetic potential A{x) such that cnT\A{x) = B{x) 
is defined up to a gauge term grad ip{x). In particular for magnetic fields (11. 5p in the class 
of potentials A{x) = (0,0, — a(r)) one can always add to a(r) an arbitrary constant c. 
This leads to the transformations Xn,m{p) ^ Xn,m{p — c) and ipn,mij,p) ^— ^ 'ipn,m{f^,P — c). 

1.3. The precise definitions of the operators and H and their decompositions into 
the direct integrals over the operators Hm{p) and H{p) are given in Section 2. To put 
it differently, we construct a complete set of eigenfunctions of the operator H. They are 
parametrized by the magnetic quantum number m, the momentum p in the direction of 
the xs-axis and the number n of an eigenvalue Xm,n{p) of the operator Hmip)- Thus, if 
we set 

i^n^ruA^.O.x,) = e^^^^e^™V„,„(r,p), (1.7) 

then 

H^„.m.p = KnAp)'^n,m,v (1-8) 

In Section 3, we show that for all ri G N and m G Z : 

• Under very general assumptions Xn,m{p) ^ oo as p ^ oo (Proposition 13.31) . 

• If 6(r) as r — > cxD, then \n,m{p) — > as p ^ — oo (Proposition 13. 5p . 

• If h{r) admits a finite positive limit as r oo, then Xn,m{p) (2^ — l)&o ioT 
all m as p —oo (Proposition 13. 6p . 

• If 6(r) ^ oo as r ^ oo, then Xn,m{p) ^ oo as p ^ — oo (Proposition 13. 6p . 

Related results concerning the dispersion curves for Schrodinger operator with con- 
stant magnetic fields defined on unbounded domains f2 C have been obtained in [5] 
(the case where is a strip) and in [4j, ^ Section 4.3] (the case where f2 is a half-plane). 

In Theorem 13.81 we formulate the main spectral results which follow from the asymp- 
totic properties of the dispersion curves Xn^m{p), p G R. First, the analyticity and the 
asymptotics as p — oo of Xn,m{p) imply immediately that the spectra a(iim) and o"(H) 
of the operators H^, m G Z, and H are purely absolutely continuous. Moreover, 

cr(H^) = [£m,oo), m G Z, o-(H) = [£o,oo), (1.9) 

where 

= inf Ai,„(p) >0. (1.10) 

Next, in the case where the magnetic field tends to as r — > oo, the spectra of H^, 
m G Z, coincide with [0, oo) and have infinite multiplicity. On the other hand, in the 
case where the magnetic field tends as r ^ oo to a positive finite limit, or to infinity. 
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we have that > for all m & Z and each of the spectra a{lim) contains infinitely 
many thresholds. 

Further, in Section 4, we obtain a convenient formula for the derivatives ^{p) which 
play the role of asymptotic group velocities. Our formula for A^^(p) yields sufficient 
conditions (see Theorem 14.31) for positivity of these functions. The leading example 
when these conditions are met, is 

b{r) = bor~', 6o>0, 5g[0,1], (1.11) 

and m 7^ 0. If 5 = 1, this result remains true for all m (cf. [I5]). On the contrary, if 5 = 
and m = 0, then A^g(p) < for all n on some interval of p (lying on the negative half- 
axis). Similar results concerning for the dispersion curves for the Schrodinger operator 
with constant magnetic field, defined on the half-plane with Dirichlet (resp., Neumann) 
boundary conditions, can be found in [1] (resp., [3] and [71 Section 4.3]). 

Finally, in Section 5 we discuss the long-time behaviour of a quantum particle. The 
time evolution of a quantum system is determined by the unitary groups exp {—iHmt), 
m G Z, so that an analysis of its asymptotics as t — > ±oo relies on spectral properties 
of the operators H^. Since these operators have discrete spectra, a quantum particle 
remains localized in the (xi, X2)-plane. Its propagation in the xs-direction is governed 
by the group velocities X'nmip)- particular, the condition X'nmip) > all n G N 
and p G M implies that a quantum particle with the magnetic quantum number m 
propagates as t ^ +oo in the positive direction of the xs-axis. 

Let us compare these results with the long-time behaviour of a classical particle in 
magnetic field (11.51) . As shown in [13], the function x'^{t) is periodic with period T 
determined by initial conditions. Its drift X3(T) — ^3(0) over the period is nonnegative 
if b{r) > and b'{r) > 0. Morever, it is strictly positive if b'{r) > for all r. In the 
case b{r) = const it is still strictly positive if the angular momentum m of a particle is 
not zero. Thus, our results for functions (II. lip correspond completely to the classical 
picture if 5 = 1 or 5 G [0, 1) and m 7^ 0. In the case 6 = and m = the behaviour of 
quantum and classical particles turn out to be qualitatively different. 

2. Hamiltonians and their diagonalizations 

Here we give precise definitions of the Hamiltonians and discuss their reductions due 
to the cylindrical symmetry. 

2.1. For an arbitrary magnetic potential A : — > such that A G L^^^^(M^)'^, the 
self-adjoint Schrodinger operator (II. 2p can be defined via its quadratic form 

h[ul=/K(Vu)W+AWuWr<fe, . = (.,,.,,.3). (2.1) 

It is easy to see that this form is closed on the set of functions u G L^(]R'^) such that 
Vu G Llc(M^f and iVu+Au G L'^iR^f. Similarly, if a G Lf^^iM."^), then the self-adjoint 
Schrodinger operator (II. 3p can be defined via its quadratic form 

h[u;p]= {\(yu){x)f + {a{x) +pf\u{x)\'^) dx, x={xi,X2), p G M. (2.2) 

JR2 
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This form is closed on the set of functions u G L^(]R^) such that integral (12.21) is finite. 
Clearly, this set does not depend on the parameter p G M. 

Let JF : L^(M^) — > L^(M^)) be the Fourier transform with respect to x^, i.e. 

{J^u){xi,X2,p) = [ e~'''''^u{xi,X2,X3)dX3. 



/27l 

If A[x) is given by formula (11. ip . then 

/oo 
h[{J='u){p)]p]dp 
-oo 

which implies the equation 

{J^B.u){xi,X2;p) = {H{p)T\i){xi,X2]p). (2.3) 

This equation can be regarded as a "working" definition of the operator H. 

2.2. Assume now that the function a in (11.11) depends only on r, and 

aGLL([0,oo);rrfr). (2.4) 

If we separate variables in the cylindrical coordinates (r, 9, x^) and denote by S^m C 
L2(M^) the subspace of functions f(r, X3)e*'"^ where f G L2(M+ x M; rrfrtixs) and m G Z 
is the magnetic quantum number, then 

L2(M=^) = 0^^. 

The subspaces S^m are invariant with respect to H so that restrictions of H on S^rn 
are related with H by formula 

H = 0H„. (2.5) 

Every can obviously be identified with the space L^(R+ x M.^rdrdx^) =: =; then 
Hm is identified with operator (11.61) . 

Quite similarly, if Tim C L2(M^) is the subspace of functions /(r)e*''"^ where / G 
L2(M_|_; rdr), then 

L2(M2) = 07^^. 

The subspaces Tim are invariant with respect to H{p) so that restrictions Hm{p) of H{p) 
on Tim are related with H{p) by formula 

H{p) = ^HM. (2.6) 

Every Ti^ can obviously be identified with the space L^(]R+; rdr) =: 7i; then Hm{p) is 
identified with operator (II. 4p . 

Let JF^ : S)m — * ^^^(IR; L^(R+; rdr)) be the restriction of JF on the subspace i^m- Then 
we have (cf. (Q) 

(.F„H„f)(r;p) = (iJ„(p) J-„f)(r; p). (2.7) 
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Sometimes it is more convenient to consider instead of Hm{p) the operator 

LUp) = r'I'HMr-"' = -f-, + + {a{r) + pf (2.8) 

acting in the space L^(M+) and unitarily equivalent to the operator Hm{p)- It is easy 
to see that the operator Lm{p) corresponds to the quadratic form 

POO 

lm[9;p] = / {\9'ir)\' + (m^ - l/A)r-'\gir)\' + (a(r) + p)2|^(r)p) rfr, (2.9) 

defined oris; inally on C^(]R+), and then closed in L2(M+). 
2.3. If 

a(r) ^ oo as r ^ oo, (2.10) 
then the spectrum of the operator Hm{p), p G M, m G Z, is discrete. Thus, it consists 
of the increasing sequence Xn,m{,p) of simple eigenvalues. Since Hm{p), p G M, is a Kato 
analytic family of type (B) (see P Chapter VII, Section 4]), all the eigenvalues Xn,m{p) 
are real analytic functions of p G M. Moreover, Xn,m{p) > because form (12.21) is strictly 
positive. 

In view of formula (12. 7p spectral analysis of the operators reduces to a study of 
a family of functions Xn,m{,p), n E N. Indeed, let An,m be the operator of multiplication 
by the function Xn,m{p) in the space L^(M). We denote by ipn,m{r;p) real normalized 
eigenfunctions (defined up to signs) of the operators Hm{p) and introduce an isometric 
mapping 

■^n,m ■■ ^^(K) L2(M+ X rdrdp) 

by the formula 
Then 



{^n,mW){p) = i)n,m{r,p)w{p). (2.11) 

L^(M+ X M; rdrdp) = Ran^„,™ 



nGN 



and 



Together with (12. 5p . formulas (12. lip and (I2.12p justify equations (II. 8p for functions (11.71) . 

3. Dispersion curves and spectral analysis 

3.1. In this subsection we consider the operators H{p) acting in the space L^(M^) by 
formula (11.30 . Under the assumption a G L^q^(]R^) they are correctly defined by their 
quadratic forms (12.21) . If 

a(x) — * oo as |x| oo, x={xi,X2), (3.1) 

then the spectrum of H{p) consists of eigenvalues A„(p), G N. We enumerate them in 
the increasing order with multiplicity taken into account. Our goal is to investigate the 
asymptotic behaviour of the eigenvalues A„(p) a.s p —>■ oo. Below we denote by C and c 
different positive constants whose precise values are of no importance. 



1/2 

v{x)\u{x)\'^dx < C sup ( / v'^{y)dy 

'\x-y\<e 
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We use the following elementary 
Lemma 3.1. Let v{x) > 0. For an arbitrary e > 0, we have the inequality 

/ / {e\Vu{x)\'^ + e-^\u{x)\^) dx 

(3.2) 

provided the supremum in the right-hand side is finite. 

Proof. Let C be a square of length e. We proceed from the estimate 

[ \u{x)\^dxY < C ( e [ \Vu{x)\^dx + e-^ f \u{x)\^dx 

which follows from the Sobolev embedding theorem by a scaling transformation. Using 
the Schwarz inequality, we deduce from this estimate that 

v{x)\u{x)\^dx < C i^j v'^{x)dx^ j \V u{x)\^ dx + j \u{x)\^dx^ . 

(3.3) 

Let us split the space in the lattice of squares Ile"^ of length e. Applying (13.31) to 
every ni"^ and summing over all n, we arrive at (13.21) . □ 

In the following assertion we do not assume (13. ip . 

Proposition 3.2. Let a G Lf^^{R^). Set a_(x) = max{— a(x), 0}, 

a{e) = sup / a^_{y)dy 

x(m? J\x-y\<e 

and suppose that ^(e) — > as e 0. Then we have 

liminfp-^inf (T(if(p)) > 1. (3.4) 

Proof. Applying estimate (13.21) with e = to the function f = a_, we find that 
/ (iVwp + (p + a)2) Iwprfx > / (iVwp + (-2pa_ +p2)|Mp) da; 

> [ (I Vm|' + \u\^dx - C^/^^ [ (I Vm|2 + p2) \u\^dx. 

Since a{p^^) ^ as p ^ oo, this implies (13.41) . □ 

Proposition 3.3. Let a G Li^^CR"^) and let condition (13.11) be satisfied. Then, for all 
G N, we have 

A„(p) =p2(i+o(l)), p^oo. (3.5) 

Proof. Under condition (13. ip the function a_ has compact support so that we can use 
Proposition 13.21 and estimate (13. 4p implies 

liminf p~^A„(p) > 1. (3.6) 
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Set G{e) = — A + {l+e~^)a^{x), e > 0. The spectrum of G{e) is discrete; let n eN, 
be the increasing sequence of its eigenvalues. By the elementary inequality 

{a + pf < {l + e-^)a^ + {l + e)p^, e > 0, 

we have H{p) < G{e) + {1 + e)p'^ so that by the minimax principle 

\n{p) < z/„(e) + (l + ey. 

Therefore, for all e > 0, 

limsupp~^A„(p) <l + e, 

which combined with (13.61) yields (13.51) . □ 

Corollary 3.4. Suppose that the function a depends on r only. Let conditions (\2A\\ 
and (I2.10p be satisfied. Then, for all n & N , m & TL, we have 

3.2. From now on we always assume that the function a depends on r only and 
that conditions (12.41) and (I2.10p are satisfied. In this subsection we investigate the 
asymptotics as p — > — oo of the eigenvalues Xn,m{p) of the operators Hm{p). Actually, 
it is more convenient to work with the operators Lm{p) acting in the space L^(]R_|_) by 
formula (12.81) . We suppose that the function a is different iable at least for suffiently big 
r and formulate the results in terms of the function b{r) = a'{r) related to the magnetic 
field by formula (II. 5p . 

Remark first that ii k = —p > is big enough, then the equation 

a(r) = k (3.7) 

has at least one solution. We denote by pk the greatest solution of (13.71) . Clearly, 

Pfc — s> oo as A; ^ oo. 

Proposition 3.5. Suppose that 

lim b{r) = 0. (3.8) 

r— >oo 

Then for each n eN and m E Z we have 

lim Xn,m{-k) = 0. (3.9) 

fc— >oo 

Proof. Set 

b(r) = sup |6(x)| and 7^ = b(pfe)-i/2. (3.10) 

x>r 

Let us fix G N. We pick a function (pi E C^(R) such that supp0i = [O, and, for 
n > 1, set 

(j)j{x) = (j)i{x - {j - l)/n), xeR, j = 2, . . . ,n. 
For k > large enough, we put 

ipj(r; k) = 7fc ^^Vj ^^^^ ^ , r>0, j = l,...,n. (3.11) 
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We will prove now that for quadratic form (12.91) 

lim lm[<^j{k); -k] = 0. (3.12) 

fc— >oo 

It follows from (13.111) that 

\^'^{r;k)\'dr<C%' (3.13) 
with C independent of k. Further, since supp (pj{k) C [pk, Pk + Ik], we have 

r~'\cp,ir;k)\'dr<Cp,\ (3.14) 

Jo 

Similarly, 

fOO 

(a(r) - kf\ipj{r; k)\'^dr < C sup (a(r) - k)^. (3.15) 
Using the condition a{pk) = k, we obtain, for r > pk, the bound 

(a(r) - kf = (a(r) - a{pk)f = b{s)ds^ < (r - pk)'h\pk) 

where b is function (l3.1Qp . Thus, the right-hand side in (13.151) is bounded by C'ylh'^{pk). 
Putting together this result with inequalities (13.131) . (I3.14p and taking into account 
(KTUlf . we get 

lm[^jiky,-k] <C{hipk)+pf) . 

This yields (Km . 

Let us use now that the supports of the functions ipj{k), j = 1, . . . , n, are disjoint and 

set 

Cn{k) = span {ipi{k), ip„{k)} . (3.16) 

Then dim Cn{k) = n and according to (I3.12p lm['^{k)] —k] — > as ^ oo for all 
^{k) G Cn{k) with ||(y9(/c)|| = 1. By the mini-max principle this implies (13.90 . □ 

The proof of Proposition 13.61 relies on a comparison of the operator Lm{—k) with the 
"model" operator 

Tik) = -^ + b\gk)ix-0kf, xeR, (3.17) 

acting in the space L^(R). Let fj be the normalized in L^(M) real- valued eigenf unctions 
(defined up to sign) of the harmonic oscillator, i.e. 

-/;(x)+a;2/,(x) = (2j-l)/,(x), xeR, jeK (3.18) 

Then 

^,{x; k) = h{QkY'^fMQkf'\^ - Qk)) (3.19) 
are normalized eigenf unctions of the operator T{k), that is 

T{k)ij,{k) = h{Qkm-l)i^Ak), JGN. (3.20) 

The proof of the following result follows the general lines of the proof of [2^, Theorem 
11.1]. 
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Proposition 3.6. Suppose that a{r) is locally semibounded from above. For r > large 
enough, we assume that the function b{r) is differentiable and that conditions 

6(r) > 0, (3.21) 
lim r'^b{r) = oo, (3.22) 

r—*oo 

as well as 

lim 6(r)-^b^(r) = 0, where bi(r) = sup \b'{x)\, (3.23) 

r/2<x<3r/2 

are satisfied. Let also 

lim k-%{pk) = 0. (3.24) 

fc— >oo 

Then, for all n & N , m & TL, we have 

Xn,m{-k) =b{Qk){2n-l + o{l)), k^oo. (3.25) 
Proof. Due to the minimax principle, it suffices to show that: 

(i) For each n G N and sufficiently large k there exists a subspace Cn{k) of L^(]R_|_) such 
that dim£„(fc) = n, Cn{k) C D{Lm{—k)), and for each ip{k) G Cn{k) we have 

{Lmi-kMk),ipik)) <biQk)i2n-l + oil))yik)f, k^oo. (3.26) 

(ii) For each n G N there exists a bounded operator Rn{k) such that rank/2„(A;) < n — 1 
(hence, Ri{k) = 0), and 

Lm{-k)>biQk)i2n-l + oil))I + Rr,{k), k^oo. (3.27) 

We pick 7fe > such that 

7fc 0, (3.28) 

Ik QkbigkY^^ ^oo, (3.29) 

7fc'&(f?fc)~'^'bi(^fc) ^0 (3.30) 

as A; ^ oo. Note that (K29h is compatible with (K2E^ due to ([Ml, and flCTj) is 
compatible with ^TMi due to flX^ . 

Proof of (i). Let C e Co°°(M) be such that < C(a;) < 1 , C{x) = 1 for |x| < 1/2 and 
suppC = 1, !]• For k large enough, set 

C{r;k) = a^kb{gky^\r-gk)), r G M+, (3.31) 

and 

ip,ir;k) = ^,{r;k)C{r;k), r G M+, j G N, (3.32) 
the functions ipj{r] k) being defined in fl3.19p . It follows from fl3.29p that 

snpp ipj{k) = [Qk - lk^b{gk)~'^^'^ , Qk + Ik^KSk)''^^'^] C [^ifc/2, 3^fc/2] 
and, in particular, fj{k) G D{Lm{—k)). Note that 

{iPj{k),i^i{k))L2(^K^) = 6ji - / ipj{x;k)4ji{x;k){l-C^{x;k))dx = 6ji + o{l) (3.33) 
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a.s k —>■ oo. Indeed, the integral here can be estimated by 



\fjix)fiix)\il - C{lkx))dx < / \f,{x)fi{x)\dx 

which tends to zero according to fl3.28p . In particular, fl3.33p implies that for all n G N 
the functions ipiik), . . . , '^n{k) are linearly independent if k is large enough. Thus, the 
space Cn{k) defined by (I3.16p has dimension n. 
Let us set 

n 

4){x; k) = '^^Cjtpjlx; k), Cj G C, (3.34) 
i=i 

ip{r; k) = ip{r; k)({r; k) and consider {Lm{—k){p{k), (p{k)). Integrating by parts, we find 
that 

-2Re {^|J\k)C{k), mem - imcik), mm) = iimcikw 

so that 

{Lmk)m,v{k)) = Re{-ij"{k) + {a{r)-kfm,^{k)C\k)) 

+ \\mc{k)r+{m'-i/mr-'mr- (3.35) 

We assume that ||v2(A;)|| = 1 and hence according to fl3.33p || = l+o(l). The second 

and third terms in the right-hand side of (13.350 are negligible. Indeed, differentiating 
(I3.3ip and using condition (13.280 . we find that 

\\mC{k)r = 0{b{Q,H) = o{b{g,)). (3.36) 

Since r^^ < 2g^^ on the support of ^p{k), relation (I3.22p implies 

\\r-'mf = Oig,') = oibigk)), fc^oo. (3.37) 

Further we consider the first term in the right-hand side of (13.350 . It follows from 
equation (I3.2O0 that 

- ^'^{k) + (a(r) - ky^,{k) = h{gu){2j - l)i;,{k) + a{k)ij,{k) (3.38) 

where the function 

a{r- k) = (a(r) - kf - h\gu){r - gtf. (3.39) 

Let us estimate the right-hand side. In view of the equation a{pk) = k, a second-order 
Taylor expansion of a at gk yields 



a{r) = k + b{gk){r - gk) + / b'{s){r - s)ds. 

Therefore, 



a{r;k) = 2b{gk){r- gk) I b'{s){r - s)ds + i I b' {s){r - s)ds' 

•> ek Qk 

and hence 

|a(r; k)\ < 6(^fe)bi(^fc)|r - ^^1^ + 4:-^hl{gk){r - guf 

< iu'KQkr'''MQk) + A-'iu%gk)-'h\{gk), 
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provided that \r — gk\ < %^b{gk)^^^'^ ■ In view of conditions fl3.28p and fl3.30p . this gives 
us the estimate 

sup \a{r;k)\=o{b{pk)) (3.40) 

|r-efe|<7,-'6(ffe)-l/' 

SO that 

{-^;{k) + (a(r) - kfij^ik)) Cik) = b{gk){2j - + o{b{p,)). 

Thus, using also fl3.33p we obtain that 

n 

Re{-r{k) + (a(r) - kf^{k),4j{k)e{k)) = b{e,) 5^(2j - l)c,ci{ip.j{k),vi{k)) + o{b{pu)) 

< b{gk)i2n-l) + oib{pk)). (3.41) 

Together with fl3.36p and (13.370 . this imphes estimate (I3.26P for each ip{k) G C-nik). 

Proof oi (ii). Let functions C e C'^(K) and r] G C°°(M) satisfy C{x) + 7f{x) = 1, 
X G M; moreover, as before, we require that < ({x) < 1 , ({x) = 1 for |x| < 1/2 and 
supp ( = [—1, 1]. By analogy with (13.31 1) set 

v{r-k) = r]{-fkb{gkY^\r- Qk)), r G R+. (3.42) 

Then we have 

('^{r;k) +r]^{r;k) = 1, r G M+. 

We proceed from the localization formula (known as the IMS formula - see e.g. [21 
Section 3.1]) 

LU-k) = ak)LU-k)ak) + r]{k)LU-kMk) - ('{kf - r]'{kf, 

where C(^)) v{k), C'{k) and ri'{k) are understood as operators of multiplication by the 
functions ({r, k), ri{r, k), ('{r, k) and //'(r, /c), respectively. According to (13.280 it follows 
from definitions (IX3B and that 

max (C'(r, kf + r/'(r, kf) = O {-flb{gk)) = o {b{gk)) , k oo. (3.43) 
Next, we check that 

r]ik)Lmi-kHk) > UkbigkWik) (3.44) 
with i/fc — ^ oo as A; oo. By virtue of the Hardy inequality 

it suffices to check that 

(a(r) - kf > Ukb{gk) (3.45) 

for 

r>gk+{2lkb{gky^^)~'=:gi^^ and r < f?, - (27,6(f?,)^/2^) -i ^. ^{-)_ (g ^g) 

According to (13.211) there exists tq such that the function a(r) is increasing for r > tq. 
Let first r > vq. Then 

|a(r) -k\ = \a{r) - a(f?fc)| > ±(a(4^^) - a(f?fc)) (3.47) 
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if ±(r — g\^'') > and r > tq. It follows from definition fl3.46p of the numbers g^i^^ that 
aiQ^t^) - aiQk) = / bis)ds = ±(27fc)-i6(f?fc)i/2 + / (bis) - h{Qk))ds. (3.48) 

J Qk J Qk 

The absolute value of the integral in the right-hand side can be estimated by 

„(±) 



ds / \h'{(j)\d(j 

Sk J Qk 



where the function bi is defined in (13.231) . By virtue of conditions (13.281) and (I3.30p 
this expression is oipf^^h^QkY/"^) as ^ oo. Therefore the absolute value of expression 
(I3.48P is bounded from below by (37fc)~^6(f)fe)^/^. Thus, for r > ro, estimate (I3.45P with 
z/fc = (37fc)^^ — i> oo is a consequence of (I3.47p . 

If r < ro, we take into account that a(r) is semibounded from above so that (a(r) — 
kY > 2~^k'^. Hence estimate (I3.47P with = 2~^k%{Qk)~^ — > oo is satisfied according 
to condition 

Putting together definitions (12. 8p and (13.170 of the operators Lm{—k) and T{k), we 
see that 

ak)L^{-k)ak) = ak)Tik)ak) + am^k), (3.49) 

where a{k) is the operator of multiplication by function ( 13.390 . The first term in the 
right-hand side is bounded from below by &(^?fc)C^(^) because h{Qk) is the first eigenvalue 
of the operator T{k). By virtue of (13.401) the second term satisfies the estimate 

\\a{k)em=o{h{Qk)). (3.50) 

It follows that operator (I3.49P is bounded from below by b{gk)C'^{k) — o{b{gk))I- Com- 
bining this result with (I3.43P and (I3.44p . we get estimate (I3.27P in the case n = 1. 

If 72 > 2, we denote by Pn{k) the orthogonal projection onto the span of the first n — 1 
eigenfunctions of the operator T{k). Then T{k){I — Pn{k)) > {2n — 1)(J — Pn{k)) and 
hence 

ak)Tik)ak) = mTmi - Pnmm + ak)nk)Pn{k)ak) 

> b{gk){2n - 1)C(A;)(/ - P„(A;))C(A;) + C(A;)T(A;)P„(A;)C(A;) = b{gk){2n - 1)C'(A;) + 

(3.51) 

where 

Rn{k) = C{k){T{k) - b{gk){2n - l)J)P„(fc)C(fc). 

Clearly, ranki?„(A;) < n - 1. Putting together ([Ml, and fl3:i9|l - fl33B . we 

obtain (I3.27P in the case n > 2. □ 

Example 3.7. Let b{r) = b^r"^ , 6o > 0, 5 < 1, for sufficiently large r. Then bi(r) = 
bo6r-^-^ and conditions (13:211) - ^TM are satisfied. Moreover, pk = cik" and k%{pk) = 
C2k^^^^ where z/ = (1 — 5)~^ and Ci,C2 > if 5 < 1. If 5 = 1, then pk = exp{bQ^k) 
and k%{pk) = fc^&o exp(— 6q ^fc). In both cases condition (13.240 is also satisfied. Thus, 
Proposition 13.61 implies the following results. If 5 > 0, then Xn,m{p) ^ as p ^ — oo 
(this result follows also from Proposition 13. 5p . If 5 = 0, then the functions Xn,m{p) have 
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finite limits 69 (2^ — 1) a.s p —00. If 5 < 0, tlien tliese functions tend to +00 as 
p —>■ —00. 

3.3. Let us return to tlie Hamiltonians and H defined in Section 2. 

Theorem 3.8. Assume ((23!) ^^^^ fl2TT0l) . 

(i) Then all operators H^, m G and hence H are absolutely continuous and their 
spectra coincide with the half-axes defined by equations fll.9p and f ll.lOp . 

(u) // the hypotheses of Proposition 1 3'. 51 hold true, then Sm = for all m G Z. 
Moreover, the multiplicities of all spectra a{H.m) and hence of a{li) are infinite. 

{Hi) Let the hypotheses of Proposition [X6l hold true. If b{r) 00, then the infimum 
in (11.101) is attained {at a finite point) so that for all m ^ Z 

£ra = minAi,„(p) > 0. 



{iv) Let the hypotheses of Proposition 1 5'. 61 hold true. If b{r) admits a finite positive 
limit 60 as r 00, then Sm G (0, 60] fo^ 0,^^ m G Z. 

Proof. It suffices to prove only the assertions concerning the operators H^. In view 
of decomposition (12.121) they reduce to corresponding statements about the operators 
A„ These operators are absolutely continuous because the eigenvalues \n,m{p) are real 
analytic functions of p G M which are non constants since according to Corollary 13.41 
^n,m{p) ^ 00 as j9 ^ 00. Moreover, we have that 

0-(An,m) = [Sn,m, Oo) whcrC £n,m = inf Xn,m{p) > (3.52) 



because Xn,m{p) > for all p G M. This implies relations (II. 9p with Sm defined by (11.101) . 

In case (ii) it suffices to use that according to (13. 9p £n,m = and hence a{An,m) = 
[0, oo) for all m and n. 

In case (iii) Proposition 13.61 implies that Xn,m{p) oo as p ^ —oo for all n and m so 
that 

Sn,m = minA„,m(p) > (3.53) 

and hence infimum in (1 1.1 OP can be replaced by minimum. 

In case (iv) we use that according to (I3.25P Sn,m < (Sri — l)6o- Moreover, £n,m > 
because Xn,m{p) > for all p G M. For n = 1, this gives the desired result. □ 

Remark 3.9. According to (I2.12p and (I3.52p the spectrum of the operator consists 
of the "branches" [Sn,m, oo) where the points £n,m are called thresholds. In cases (iii) 
and (iv) 

£n,m ^ £n+l,m (3.54) 

for all n G N. Indeed, in case iii) (I3.54p is a consequence of the estimate Xn^m{p) < 
Xn+i,m{p) valid for all p G M and of formula (I3.53p . In case (iv) one has to take addi- 
tionally into account that the limit of Xn,m{p) as p ^ — oo is strictly smaller than that 
of Xn+i,m{p)- Inequality (I3.54p means that there are infinitely many distinct thresholds 
in each of the spectra a{H.m)i m E and hence in o"(H). 



TRANSLATIONALLY INVARIANT MAGNETIC OPERATORS 



15 



Remark 3.10. In case( iii) the multiplicity of the spectrum of all operators A„ m equals 
at least to 2 whereas in cases (ii) and (iv) it might be equal to 1. 

4. Group velocities 

4.1. In this subsection we obtain a formula for the derivative ^(p), n G N, m G Z, 
which yields sufficient conditions for the monotonicity of \n,m{p) as a function of p. 
Recall that the operators Hm{p), m E N, p E M, were defined in the space H by formula 

The proof of Theorem 14.31 relies on integration by parts. To prove that non-integral 
terms disappear at r = 0, we use standard bounds on ipn,m{r; p). Unfortunately, we were 
unable to find necessary results in the literature and therefore give their brief proofs. 

Let us consider the differential equation of Bessel type 

-r'^{ry'y + m^r''^y + q{r)y = m = 0,1,2,..., (4.1) 

in a neighborhood (0,ro) of the point r = 0. If g(r) = 0, then it has the regular 
yQ^^\r) = r"^ and singular y^^^^\r) = r"*^ solutions for m ^ and yQ^^\r) = 1 and 
y^^^^\r) = Inr for m = 0. 

Lemma 4.1. Let m ^ 0, and let the function rq{r) belong to the class L^(0,ro). Then 
equation (14. ip has a solution y(^^^)(r) satisfying the relation 

?/("'^»)(r) = r™ + o(r™), r ^ 0. (4.2) 

For its derivative, we have the bound 

dy^'''^\r)/dr = Oir"""^). (4.3) 

Let m = 0. Suppose that the function rlnrg(r) belongs to the class L^(0,ro). Then 
equation (14. ip has a solution ?/(''"^5)(r) satisfying relation (14.20 where m = 0. For its 
derivative, we have the bound 

dy^''^^\r)/dr = o(^j^ \q{s)\d^ . (4.4) 

Moreover, if the function r\v?rq{r) belongs to the class L^(0,ro), then equation (14. ip 
has a solution y(^*"5)(r) satisfying the relation 

?/("*"f)(r) = Inr + o(l), r ^ 0. (4.5) 

In this case any bounded solution of equation (14.10 coincides [up to a constant factor) 
with the regular solution y^^'^^\r). 

Proof. We construct the function y(^^5)(r) as the solution of the Volterra integral equa- 
tion 

y^^^'\r) = y;[^'\r) + K^ F s{yt'\r)yl,-'\s)-yl;^'\s)yi-'\r (4.6) 

^0 
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where = (2m) ~^ for m ^ and xq = —1. Differentiating it twicely, we see that 
yireg)^^-^ satisfies equation (14 .11) . Equation (14.61) can be solved by iterations, that is 

oo 

y^'''\r) = J2yl:''\r). (4.7) 

n=0 

Hereby the n^'^-iteration obeys the bound 

\y^:''\r)\<—r^( / s\qis)\ds 
^- ^ Jo 

if m 7^ 0; if m = 0, then s|g(s)| should be replaced by s\ lns||g(s)|. This ensures the 
convergence of series (14. 7p as well as relation (14. 2p . Differentiating equation (14. 6p and 
using (14. 2p . we get bounds (14. 3 p and (14.40 on the derivative of y^^^^\r). 

If m = 0, we can construct the function ?/(**"9)(r) as the solution of equation (14. 6 p 
where the first term, y^'^^\r), in the right-hand side is replaced by y\^''^^\r)., that is 

rr 

(«^«9)(r) = Inr + / s\n{r / s)q{s)y^''''^\s)ds. 



This equation can again be solved by iterations which, in particular, implies estimate 

m. □ 



This result can be supplemented by the following 

Lemma 4.2. Let m ^ 0, and let the function rq^{r) belong to the class L^(0,ro). 
Assume additionally that q = q. If ip is a solution of equation (14.11) from the class 
L^((0, To); rrfr), then it coincides {up to a constant factor) with the regular solution 
y^^'^9)[r) and hence satisfies estimates (14.21) and ( 14. 3p . 



Proof. Let us extend the function q{r) to (ro, oo) by zero, and let us consider the differ- 
ential operator 

hy = —r^^iry')' + nn?r~'^y + q{r)y 

in the space L^(]R+;r(ir) on domain C^(M+). If g = 0, we denote this operator by h^. 
The operator Kq is essentially self-adjoint. To prove the same for h, it suffices to check 
that 

/ q\r)\f{T)\\dT<e\\h,fr + C\\f\\\ / G Co~(M+), e<l. (4.8) 
Let us use the estimate 

q^{\x\)\u{x)\^dx < I g^(|x|)(ix max |M(a;)|^ 

x\<ro J\x\<n) '-^^"^"^ 

<e [ \{Au){x)\^dx + Ce~^ [ \u{x)\^dx, We > 0. 

Restricting it on the subspace of functions u{x) = /(r)e*™^, we obtain estimate (14. 8 p 
which implies that h is essentially self-adjoint as well as h^. Thus, equation (14. ip has at 
most one solution from L'^{{0, Tq); rdr) which is necessarily proportional to y^^'^^\r). □ 
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Now we are in a position to obtain a formula for the derivative A^^(p). In addition 
to our usual assumptions that b{r) is not too singular at r = 0, an integration-by- 
parts marchinery requires that b{r) does not vanish too rapidly as r ^ 0. The precise 
conditions are formulated rather differently in the cases m ^ and m = 0. We start 
with the first case. 

Theorem 4.3. Let m ^ 0. Suppose that b G C^(M_|_) and b{r) > 0, r E IR+. Assume 
fl2.10p and that b{r) = 0{e'^^) for some c>Oasr— >oo. At r = we suppose that 
b{r) = 0{r~"') where 7 < 3/2. Moreover, we assume that for some (3 < 2\m\ — 1 

< A; = 0,1, 2, 3, (4.9) 

Put 

v{r) = r{r~^{rb{r)~^yy . 

Then 

PCO 

A:,„(p) = -2 / rb~'ir)b'irU^^^ir;pfdr 
Jo 

POO poo 

-2-1 / v'ir)ijl^ir;p)dr + 2m'' / r-%-\r)^l^ir;p)dr, (4.10) 
Jo Jo 

where the eigenfunctions ipn,m{f'',p) of the operator Hm{p) are real and normalized, that 

is \\tpn,m\\ = 1- 

Proof. In view of the equation 

(a(r) + pfi)n,m = {r^lj'^,m)' - m^r~^i^n,m + K,mi^n,m (4.11) 

we can apply to the function ipn,m the results of Lemmas 14.11 and 14.21 where q{r) = 
(a(r) +py — \n,m- Thus, Lemma implies that '4'n,m{r'-iP) = 0(r'™l) and '^'nmi'^'^P) — 
0(rl™l~i) as r — s> which ensures that non-integral terms disappear at r = 0. 

To prove the same for non-integral terms corresponding to r oo, we use super- 
exponential decay of eigenfunctions ipn,m{r; p) of the operators Hm{p)- This result is valid 
[H] (see also [6J) for all one-dimensional Schrodinger operators with discrete spectra. In 
view of the condition a(r) = 0{e'^^), it follows from equation (14.111) that the derivatives 
i^nmi^'^P) ^Iso decay super-exponentially. 

Let us proceed from the formula of the first order perturbation theory (known as the 
Feynman-Hellman formula) 

rd(a(r)+p)\o , , , P 5(a(r-) +p)2 ,2 . , , ^ 
K,m{p) = l ^l^{r;p)rdr = -±±l-J^^l^{r;p)T{r)dr (4.12) 

where r(r) = rb{r)~^. Using that a(r) = 0{r^^"') and r(r) = 0{r^~^), we integrate by 
parts and get 

POO 

K,mip) = - {a{r) + pYtpn,m{r;p){T'{r)tpri,m{r;p) + 2T{r)tp'^„^{r;p))dr. 
Jo 
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Now it follows from equation (14.111) that 

K,mip) = -Kmip) I {'r{r)^l^^{r]p))'dr + w? I r^^{T{r)%l)l^^{r]p))' dr 
Jo 

r'^i^i)' {r; p)y{T'{r)^n,m{r; p) + 2r(r)V^' (r; p))dr. (4.13) 



By the condition T{r)ilj^^{r;p) — as r — 0, the first term in the right-hand side 
equals zero. In the second term we integrate by parts which yields 

POO 

r-'^{T{r)^l {r;p)ydr = 2 / r-^T{r)iljl (r;p)dr 



because r~^r(r)?/^^ ^(r; p) — > 0. 

In the last integral in the right-hand side of (I4.13p . we also integrate by parts using 
that i'n^mi^j P)i^n,m{r; p)t' (r) ^ as r — > 0. Thus, we have that 

oo POO 



r {ri^n,mir'^P))'^'{r)4jn,m{r;p)dr = / r'(r)^;^(r;p) dr 

Jo 

POO 

+ / v{r)ip'^^^{r;p)ipn,mir;p)dr(4:.U) 
Jo 

The last integral in the right-hand side equals 

POO 

-2-^ / v'{r)iljl^{r;p)dr 
Jo 

because v{r)ilj^^{r;p) ^ as r — >• 0. Similarly, we get that 

POO POO 

-2 / r-\rilj' (r■,p)yT{r)^p' {r;p)dr = - r-^T{r)d{rip' (r;p)'^) 



r (r r(r))V;,„(r;p) dr 

since T{r)ip'^^{r;pY — > as r — 0. Putting the results obtained together, we arrive at 
representation (I4.10p . □ 

Corollary 4.4. Ifh'{r) < and r%{r)v'{r) < Am? for all r > 0, then A^^(p) > for 
allp&M. and n. If, moreover, one of these inequalities is strict on some interval, then 

K,miP) > 0- 

Corollary 4.5. If b{r) = bor'^ 6 e [0, 1], then r(r) = bQ^r^+\ v{r) = b^^S"^ - ly-^ 
and 



A'„,^(p) = 26o ^5 / rX„(r;p)2rfr 

^0 

POO 

+b^\2m' - 2-1(1 - Sfil + 6)) / (r;p)rfr. 
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For bo > 0, this expression is strictly positive {so that the functions Xn,m{p) ore strictly 
increasing for all p G M) for m ^ since (1 — 5)^(1 + 5) <1. Moreover, for 6 = 1 this 
result is true for all m E 1j. 

In the case m = we consider for simplicity only fields (II .lip . 

Proposition 4.6. If b{r) = h^r^^ , 5 E [0, 1], then 

POO 

A;,o(p) = 2&0 ' W r'^P'^,oir;prdr 
Jo 

POO 

-b,'2-\l - + 6) / r-'^\4^lo{r;p) -^loiO; P))dr. 



Ifbo>0 and 6 = 1, then A^ o(p) > for all p e M. 

Proof. Let us proceed again from formula fl4.12p . We use now that the function ipnfiilxl', p) 
of X G belongs to the Sobolev class Hf^^iM."^), and therefore ipn,o{r;p) has a finite limit 
as r — s> 0. Thus, by Lemma l4.ll ■t/'^ gfr: p) = ©(r^""^) for any e > as r — 0. These 
results allow us to intergrate by parts as in the case m ^ 0. The only difference is with 
the second integral in the right-hand side of fl4.14p . Now v{r) = bQ^{5'^ — l)r^~^ and 
this integral equals 



v{r)ip'n,o{r;p)'^n,o{r;p)dr = 2W f (r)rf(V^^_o('^; P) " ^n,o(0; P)) 

Jo 

POO 

= -2-W t;'(r)«o(r;p)-<o(0;p))dr 
Jo 

because v{r)(ilj'^ Q(r; p) — ip'^ q(0; p)) as r — > 0. □ 

4.2. In this subsection we show that for linear potentials, that is for magnetic fields 
not depending on r, all eigenvalues Xnfi{p), n G N, of the operator Ho{p) are not monot- 
onous functions of p G M. We follow closely the proof of the first part of Proposition [321 
However we now use that eigenf unctions of the harmonic oscillator decay faster than 
any power of r~^ at infinity (actually, they decay super-exponentially). 

Proposition 4.7. Assume that for sufficiently large r 

b{r) = bo>0. (4.15) 

Then, for all n eN, some 7„ > and sufficiently large k > 0, we have 

Xn,oi-k) < {2n - l)6o - lnk~^. (4.16) 

Proof. Let be the same function as in the proof of the first part of Proposition 13.61 We 
set pk = b^^k, 7fc = 2b]l'^k~^ and define the functions C,{r] k) and ^j{r; k) by formulas 
fl3.3ip and fl3.32p . respectively. It suffices to check that 

(Lo(-A;)v9(A;), ^{k)) <2n~l- 7„fc-2. (4.17) 

for sufficiently large k and all normalized functions from subspace (13.161) . Let us proceed 
from formula fl3.35p . Since the functions ipj{x; k) decay faster than any power of \x\~^ 
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as |x| oo, the term o(l) in fl3.33p is actually 0{k °°). Similarly, estimate (13.361) can 
be formulated in a more precise form as 

\\mc{k)r = oik-n. (4.18) 

Since r < 2^^3k on the support of ip{k), we have that 

\\r-'tlj{k)f > {2/3yk-^. (4.19) 

Now function (I3.34p is zero if r and k are large enough. Therefore equation (I3.38P yields 
the exact equality 

n 

Re{-r{k) + {bor - k)'^{k),^{k)C{k)) = 6o $^(2j - l)cMfjik),Mk)) 

(cf. (13.411) ). Up to terms 0{k~'^), the right-hand side here is estimated by 6o(2^ ~ !)• 
Together with (1418]) and (147191) . this implies estimate KTT^ . □ 

Combining relations (13.250 and (14.160 . we see that the eigenvalues Xn,o{.p) tend as 
p — oo to their limits (2n — l)6o from below. On the other hand, according to (13. 5p 
^nfiip) ^ cxD as p — > oo. Thus, all functions A„^o(p) have necessarily local minima. We 
can obtain an additional information using the following elementary 

Lemma 4.8. Suppose that (14.151) is satisfied for all r > and that a{r) = bor. Then 

A„,™(0) =26o(2n-l + |m|) (4.20) 

for all n and m G Z. 

Proof. Let us consider the two-dimensional harmonic oscillator T = —A + b^xf + x^). 
Separating the variables xi, X2, we see that its spectrum consists of the eigenvalues 
26o(^i + ^2 — 1) where /i,/2 € It follows that the operator T has the eigenvalues 
2^0 J) j ^ of multiplicity j. On the other hand, separating the variables in the polar 
coordinates, we see that the spectrum of T consists of the eigenvalues A„_m(0) of the 
operators Hm{0)- For the proof of (14.201) we take into account that all eigenvalues 
An,m(0) are simple and that A„_m+i(0) > A„^m(0) for all n and m > 0. Clearly, the 
operator -^0(0) has an eigenvalue 26oi if cind only if its multiplicity j is odd. This gives 
formula (I4.20p for m = 0. We shall show that for every j G N 

Ai,,-i(0) = A2,,_3(0) = . . . = A2,-,+3(0) = Ai,_,+i(0) = 2boj (4.21) 

which is equivalent to formula (14.200 for all m. Let us choose some jo and suppose that 
(I4.2ip holds for all j < Jq. Then we check it for j = jo + 1. First we remark that if an 
operator Hm{0) for some m > has n eigenvalues in the interval [2bo, 26o(jo + 1)]; then 
the operator Hm-i{0) has at least n eigenvalues in the interval [26o,26ojo]- Then using 
(I4.2ip for j < jo, we see that if an operator Hm{0) has the eigenvalue 26o(jo + 1)) then 
necessarily the operator Hm-i{0) has the eigenvalue 26ojo- Therefore according to (I4.2ip 
for j = jo, only the operators i^m(O) with m = jo, jo — 2, ... , —jo -|- 2, —jo might have 
the eigenvalue 26o(jo + !)• There are jo + 1 of such operators and the multiplicity of this 
eigenvalue equals jo + 1. Thus, all the operators Hm{0) for m = jo, jo — 2, . . . , — jo + 2, —jo 
and only for such m have the eigenvalue 26o (jo + 1 ) • This proves (14.210 for jo + 1 . □ 
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Comparing this result with fl3.27p . we see that, for potentials a{r) = h^r, 
lim Xn,m{p) = &o(2n - 1) < 26o(2n - 1 + |m|) = A„.m(0). 

Together with fl4.16p . this implies that the functions Xn,o{p) have negative local minima. 
Thus, we get the following 

Theorem 4.9. Under the hypotheses of Proposition \4^ the eigenvalues A„,o(p)j n ^N, 
of the operator Ho{p) are not monotonous functions o/p G M. Moreover, if fl4.15p is 
satisfied for all r > and a{r) = h^r, then the functions \nfl{p) lose their monotonicity 
for p < 0. 

We do not know how many minima have the functions A„^o(p)- 

The problem of monotonicity of the eigenvalues A„^o(p) for fields h{r) = h^r^^ where 
5 E (0, 1) remains also open. 

4.3. In a somewhat similar situation the break down of monotonicity of group 
velocities was exhibited in [7|. In this paper one considers the Schrodinger opera- 
tor H*^^) = — ^ + ~ ^-^j with constant magnetic field 6 > 0, defined on the 

semi-plane {(x, ?/) G : x> 0} with the Neumann boundary condition at x = 0. Let 
H^^\p) = —d^/dx^ + {hx +pY, p G M, be the self-adjoint operator in the space L^(M+) 
corresponding to the boundary condition m'(0) = 0. Then the operator H^^^ is unitarily 
equivalent under the partial Fourier transform with respect to y, to the direct integral 
H^^\p)dp. It is shown in [7, Section 4.3] that the lowest eigenvalue of H^'^\p) 
is not monotonous for p < 0. This follows from the inequality /i'i(O) > provei£| in [3] 
and the relations 

lim = /ii(0) = b. (4.22) 

Our proof of non- monotonicity of the functions Xn,o{p) is essentially different since in 
contrast with fl4.22p we have limp_^_oo A„^o(p) < A„,o(0). 

5. Asymptotic time evolution 

5.1. Combined with the stationary phase method, the spectral analysis of the op- 
erators H = H(a) allows us to find the asymptotics for large t of solutions u{t) = 
exp(— iHt)'Uo of the time dependent Schrodinger equation. It follows from (11.20 that 

exp(— zH(a)t)Mo = exp(zH(— a)t)lZo. 

Therefore it suffices to consider the case a(r) — > +oo. Moreover, on every subspace 9)rn 
with a fixed magnetic quantum number m, the problem reduces to the asjTiiptotics of 
the function u(t) = exp{—i'H.mt)uo. 

Let us proceed from decomposition (I2.12p . Suppose that J-'mUo G Ran\E^„ ^- Then 
(see (EH])) 

iJ^niUo)ir,p) = i>n,mir,p)f{p) (5.1) 



Note that in [31 and [7j the parameter p is chosen with the opposite sign. 



22 



D. YAFAEV 



where / = K,m^mUo and u{t) = ^4^„,™e-^^"^-*/, that is 

/oo 
e*""^-^^"-(^)V„,„(r,p)/(p)f/p. (5.2) 
-co 

The analytic function X'n^^ip) might have only a countable set of zeros Pn,m,i with 
possible accumulations at ±00 only. The function X'nmip) monotone on every in- 
terval {pn,m,i,Pn,m,i+i) and takcs there all values between A^.^(p„,m,«) =: C(n,m,i and 
X'n miPn,m,i+i) ='■ Pn,m,i- We cousider the asymptotics of integral (15.21) on each of the 
subspaces L'^{pn,m,hPn,m,i+i) Separately. Let us set 7 = x^t'^. First we suppose that 
/ G {pn,m,h Pn,m,i+i) ■ The Stationary points of integral (15.21) are determined by the 
equation 

If 7 ^ {an,m,i,f3n,m,i), it does not have solutions from the interval {pn,m,i,Pn,m,i+i)- 
Therefore integrating directly by parts, we find that function (15.21) decays in this region 
of Xs/t faster than any power of (l^sl + Itl)""*^ (and r). If 7 G {an,m,h Pn,m,i), then on 
the interval {pn,m,h Pn,m,i+i) equation (15.31) has a unique solution which we denote by 
i^n,m,i{l)- Let us set 

^n,m,l{l) = ^n,mAl)l " (z/„,m,i (7) ) 

and denote by Xn,m,i the characteristic function of the interval {an,m,h Pn,m,i)- For 7 from 
this interval, we apply the stationary phase method to integral (15. 2p which yields 

^(r,X3,t) = r(2,e^*"-'WVn,™(r,z.„,^,K7))|A:,^(z/„,^,^ 

xfi^n,m,l{l))Xn,m,l{'y)\''^\'^^'^ + ^ooir,X3,t), 7 = Xgt'S t ±00, (5.4) 

where r^^ ^ = e=F'^*'^g'^(^".'-(p))/'^ for p G {Pn,m,h Pn,m,i+i) and 

lim \\uoo{;t)\\=0. (5.5) 

t— »±oo 

Since the norm in the space of the first term in the right-hand side of (15. 4p equals 
the norm of / in the space L"^ {pn,m,i, Pn,m,i+i) , asymptotics (15.41) extends to all functions 
(15. ip with an arbitary / G L^{Pn,m,hPn,m,i+i)- Thus, we have proven 

Theorem 5.1. Assume (12.40 and (I2.10p . Let u{t) = exp{—i'H.mt)uQ where Uq satisfies 
(I5.jp with f G L'^{Pn,m,i, Pn,m,i+i) ■ Then the asymptotics as t ^ ±oo of this function is 
given by relations ( |5.^p , (I5.5p . 

Of course asymptotics (15. 4p . (15. 5p extends automatically to all / G L^(M) with com- 
pact support and to linear of functions 4'n,m{r,p)fn{p) over different n. 

By virtue of formulas (15.40 . (15. 5p a quantum particle in magnetic field (11.50 remains lo- 
calized in the (xi, a:2)-plane but propagates in the Xs-direction. If / G L'^{pn,m,i, Pn,m,i+i), 
then a particle "lives" as \t\ — >• cxo in the region where x^ G {c(n,m,it, Pn,m,it)- In par- 
ticular, if X'{p) > {X'{p) < 0) for p G {pn,m,i,Pn,m,i+i), then a particle propagates 
in the positive (negative) direction as t +oo. Thus, according to Corollary 14.51 if 
6(r) = b^r''^, 6 G [0,1], bo > 0, then a particle with the magnetic quantum number 
m 7^ propagates always in the positive direction of the xa-axis. If 5 = 1, then this 
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result remains true from all m. On the contrary, if 5 = and m = 0, then a particle 
will propagate in a negative direction for some interval of momenta p. 

5.2. Theorem 15.11 implies the existence of asymptotic velocity in the X3-direction. 
The corresponding operator is defined by the equation (cf. fl2.12p ) 

where ^ are the operators of multiplication by the functions A'„ ^(p). To put it differ- 
ently, the operator acts as multiplication by A^^(p) in the spectral representation 
of the operator where it acts as multiplication by the functions Xn,m{p)- 

Proposition 5.2. Assume (12.41) and (12.101) . Then, for an arbitrary bounded function 

Q, 

s-lim|i|^oo exp (iH„t) Q{x^/ 1) exp (-iH^t) = Q (H'^) (5.6) 
{in particular, the strong limit in the left-hand side exists). 

Proof. We shall check that for all Uq G Tim 

lim ||Q(x3/t) exp (-iHmt)Mo - exp (-iHmt)Q(H^)Mo|| = (5.7) 

which is equivalent to relation (15.61) . Remark that if uq satisfies (15.11) . then 

{^n,mQ{^JUo){r,p) = ^PrrAr,P)Q{K,m{p))f{p)- (5-8) 

It suffices to prove (15.71) on a dense set of elements Uq such that equality (15.11) is true 
with / G L'^{Pn,m,i,Pn,m,i+i)- Applying the operator Qi^Xs/t) to asymptotic relation 
(15. 4p . we see that the asymptotics of Q^x^/t) exp {—i'Hmt)uo is given again by formula 
(15. 4p where the function f {yn,m.,i{,l)) in the right-hand side is replaced by the func- 
tion Q{'^) f {i'n,m,i{l)) ■ Similarly, it follows from Theorem 15.11 and relation (15.80 that 
the asymptotics of exp (— 'iHmt)Q(H^)-uo is given by formula (15. 4p where the function 
f{^n,m,i{l)) in the right-hand side is replaced by the function Q{\n^rn{^n,m,iil))) f ij^n,m,i{^)) ■ 
So for the proof of (15.71) . it remains to take equation (15.31) into account. □ 

Relation (15.61) shows that can naturally be interpreted as the operator of asymp- 
totic velocity in the Xa-direction. 

Similar results concerning the Iwatsuka model (see [8] or [2]) have been obtained in 

m- 

Numerous useful discussions with Georgi Raikov as well as a financial support by the 
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